Depletion of TBC1D4 improves the metabolic exercise response by overcoming genetically induced peripheral insulin resistance.

The RabGTPase-activating protein (RabGAP) TBC1D4 (=AS160) represents a key component in the regulation of glucose transport into skeletal muscle and white adipose tissue (WAT) and is therefore crucial during the development of insulin resistance and type-2 diabetes. Increased daily activity has been shown to be associated with improved postprandial hyperglycemia in allele carriers of a loss-of-function variant in the human TBC1D4 gene. Using conventional Tbc1d4-deficient mice (D4KO) fed a high-fat diet (HFD), we show that already a moderate endurance exercise training leads to substantially improved glucose and insulin tolerance and enhanced expression levels of markers for mitochondrial activity and browning in WAT from D4KO animals. Importantly, in vivo and ex vivo analyses of glucose uptake revealed increased glucose clearance in interscapular brown adipose tissue (iBAT) and WAT from trained D4KO mice. Thus, chronic exercise is able to overcome the genetically induced insulin resistance caused by the Tbc1d4-depletion. Gene variants in TBC1D4 may be relevant in future precision medicine as determinants of exercise response.

Proteomic and transcriptomic profiling of brainstem, cerebellum and olfactory tissues in early- and late-phase COVID-19.

Neurological symptoms, including cognitive impairment and fatigue, can occur in both the acute infection phase of coronavirus disease 2019 (COVID-19) and at later stages, yet the mechanisms that contribute to this remain unclear. Here we profiled single-nucleus transcriptomes and proteomes of brainstem tissue from deceased individuals at various stages of COVID-19. We detected an inflammatory type I interferon response in acute COVID-19 cases, which resolves in the late disease phase. Integrating single-nucleus RNA sequencing and spatial transcriptomics, we could localize two patterns of reaction to severe systemic inflammation, one neuronal with a direct focus on cranial nerve nuclei and a separate diffuse pattern affecting the whole brainstem. The latter reflects a bystander effect of the respiratory infection that spreads throughout the vascular unit and alters the transcriptional state of mainly oligodendrocytes, microglia and astrocytes, while alterations of the brainstem nuclei could reflect the connection of the immune system and the central nervous system via, for example, the vagus nerve. Our results indicate that even without persistence of severe acute respiratory syndrome coronavirus 2 in the central nervous system, local immune reactions are prevailing, potentially causing functional disturbances that contribute to neurological complications of COVID-19.

Neurofibromin 1 controls metabolic balance and Notch-dependent quiescence of murine juvenile myogenic progenitors.

Patients affected by neurofibromatosis type 1 (NF1) frequently show muscle weakness with unknown etiology. Here we show that, in mice, Neurofibromin 1 (Nf1) is not required in muscle fibers, but specifically in early postnatal myogenic progenitors (MPs), where Nf1 loss led to cell cycle exit and differentiation blockade, depleting the MP pool resulting in reduced myonuclear accretion as well as reduced muscle stem cell numbers. This was caused by precocious induction of stem cell quiescence coupled to metabolic reprogramming of MPs impinging on glycolytic shutdown, which was conserved in muscle fibers. We show that a Mek/Erk/NOS pathway hypersensitizes Nf1-deficient MPs to Notch signaling, consequently, early postnatal Notch pathway inhibition ameliorated premature quiescence, metabolic reprogramming and muscle growth. This reveals an unexpected role of Ras/Mek/Erk signaling supporting postnatal MP quiescence in concert with Notch signaling, which is controlled by Nf1 safeguarding coordinated muscle growth and muscle stem cell pool establishment. Furthermore, our data suggest transmission of metabolic reprogramming across cellular differentiation, affecting fiber metabolism and function in NF1.

Single-cell transcriptomic profiling reveals a pathogenic role of cytotoxic CD4+ T cells in giant cell arteritis.

Giant cell arteritis (GCA) is a systemic vasculitis mediated by an aberrant immunological response against the blood vessel wall. Although the pathogenic mechanisms that drive GCA have not yet been elucidated, there is strong evidence that CD4+ T cells are key drivers of the inflammatory process occurring in this vasculitis. The aim of this study was to further delineate the role of CD4+ T cells in GCA by applying single-cell RNA sequencing and T cell receptor (TCR) repertoire profiling to 114.799 circulating CD4+ T cells from eight GCA patients in two different clinical states, active and in remission, and eight healthy controls. Our results revealed an expansion of cytotoxic CD4+ T lymphocytes (CTLs) in active GCA patients, which expressed higher levels of cytotoxic and chemotactic genes when compared to patients in remission and controls. Accordingly, differentially expressed genes in CTLs of active patients were enriched in pathways related to granzyme-mediated apoptosis, inflammation, and the recruitment of different immune cells, suggesting a role of this cell type in the inflammatory and vascular remodelling processes occurring in GCA. CTLs also exhibited a higher clonal expansion in active patients with respect to those in remission. Drug repurposing analysis prioritized maraviroc, which targeted CTLs, as potentially repositionable for this vasculitis. In addition, effector regulatory T cells (Tregs) were decreased in GCA and showed lower expression of genes involved in their suppressive activity. These findings provide further insights into the pathogenic role of CD4+ T cells in GCA and suggest targeting CTLs as a potential therapeutic option.

Non-classical circulating monocytes expressing high levels of microsomal prostaglandin E2 synthase-1 tag an aberrant IFN-response in systemic sclerosis.

Systemic sclerosis (SSc) is a complex disease that affects the connective tissue, causing fibrosis. SSc patients show altered immune cell composition and activation in the peripheral blood (PB). PB monocytes (Mos) are recruited into tissues where they differentiate into macrophages, which are directly involved in fibrosis. To understand the role of CD14+ PB Mos in SSc, a single-cell transcriptome analysis (scRNA-seq) was conducted on 8 SSc patients and 8 controls. Using unsupervised clustering methods, CD14+ cells were assigned to 11 clusters, which added granularity to the known monocyte subsets: classical (cMos), intermediate (iMos) and non-classical Mos (ncMos) or type 2 dendritic cells. NcMos were significantly overrepresented in SSc patients and showed an active IFN-signature and increased expression levels of PTGES, in addition to monocyte motility and adhesion markers. We identified a SSc-related cluster of IRF7+ STAT1+ iMos with an aberrant IFN-response. Finally, a depletion of M2 polarised cMos in SSc was observed. Our results highlighted the potential of PB Mos as biomarkers for SSc and provided new possibilities for putative drug targets for modulating the innate immune response in SSc.

IsoTools: a flexible workflow for long-read transcriptome sequencing analysis.

MOTIVATION: Long-read transcriptome sequencing (LRTS) has the potential to enhance our understanding of alternative splicing and the complexity of this process requires the use of versatile computational tools, with the ability to accommodate various stages of the workflow with maximum flexibility. RESULTS: We introduce IsoTools, a Python-based LRTS analysis framework that offers a wide range of functionality for transcriptome reconstruction and quantification of transcripts. Furthermore, we integrate a graph-based method for identifying alternative splicing events and a statistical approach based on the beta-binomial distribution for detecting differential events. To demonstrate the effectiveness of our methods, we applied IsoTools to PacBio LRTS data of human hepatocytes treated with the histone deacetylase inhibitor valproic acid. Our results indicate that LRTS can provide valuable insights into alternative splicing, particularly in terms of complex and differential splicing patterns, in comparison to short-read RNA-seq. AVAILABILITY AND IMPLEMENTATION: IsoTools is available on GitHub and PyPI, and its documentation, including tutorials, CLI, and API references, can be found at https://isotools.readthedocs.io/.

BMP4 triggers regulatory circuits specifying the cardiac mesoderm lineage.

Cardiac lineage specification in the mouse is controlled by TGFß and WNT signaling. From fly to fish, BMP has been identified as an indispensable heart inducer. A detailed analysis of the role of Bmp4 and its effectors Smad1/5, however, was still missing. We show that Bmp4 induces cardiac mesoderm formation in murine embryonic stem cells in vitro. Bmp4 first activates Wnt3 and upregulates Nodal. pSmad1/5 and the WNT effector Tcf3 form a complex, and together with pSmad2/3 activate mesoderm enhancers and Eomes. They then cooperate with Eomes to consolidate the expression of many mesoderm factors, including T. Eomes and T form a positive- feedback loop and open additional enhancers regulating early mesoderm genes, including the transcription factor Mesp1, establishing the cardiac mesoderm lineage. In parallel, the neural fate is suppressed. Our data confirm the pivotal role of Bmp4 in cardiac mesoderm formation in the mouse. We describe in detail the consecutive and cooperative actions of three signaling pathways, BMP, WNT and Nodal, and their effector transcription factors, during cardiac mesoderm specification.

Odd skipped-related 1 controls the pro-regenerative response of fibro-adipogenic progenitors.

Skeletal muscle regeneration requires the coordinated interplay of diverse tissue-resident- and infiltrating cells. Fibro-adipogenic progenitors (FAPs) are an interstitial cell population that provides a beneficial microenvironment for muscle stem cells (MuSCs) during muscle regeneration. Here we show that the transcription factor Osr1 is essential for FAPs to communicate with MuSCs and infiltrating macrophages, thus coordinating muscle regeneration. Conditional inactivation of Osr1 impaired muscle regeneration with reduced myofiber growth and formation of excessive fibrotic tissue with reduced stiffness. Osr1-deficient FAPs acquired a fibrogenic identity with altered matrix secretion and cytokine expression resulting in impaired MuSC viability, expansion and differentiation. Immune cell profiling suggested a novel role for Osr1-FAPs in macrophage polarization. In vitro analysis suggested that increased TGFß signaling and altered matrix deposition by Osr1-deficient FAPs actively suppressed regenerative myogenesis. In conclusion, we show that Osr1 is central to FAP function orchestrating key regenerative events such as inflammation, matrix secretion and myogenesis.

Combinatorial effects on gene expression at the Lbx1/Fgf8 locus resolve split-hand/foot malformation type 3.

Split-Hand/Foot Malformation type 3 (SHFM3) is a congenital limb malformation associated with tandem duplications at the LBX1/FGF8 locus. Yet, the disease patho-mechanism remains unsolved. Here we investigate the functional consequences of SHFM3-associated rearrangements on chromatin conformation and gene expression in vivo in transgenic mice. We show that the Lbx1/Fgf8 locus consists of two separate, but interacting, regulatory domains. Re-engineering of a SHFM3-associated duplication and a newly reported inversion in mice results in restructuring of the chromatin architecture. This leads to ectopic activation of the Lbx1 and Btrc genes in the apical ectodermal ridge (AER) in an Fgf8-like pattern induced by AER-specific enhancers of Fgf8. We provide evidence that the SHFM3 phenotype is the result of a combinatorial effect on gene misexpression in the developing limb. Our results reveal insights into the molecular mechanism underlying SHFM3 and provide conceptual framework for how genomic rearrangements can cause gene misexpression and disease.

Primary ChAdOx1 vaccination does not reactivate pre-existing, cross-reactive immunity.

Currently available COVID-19 vaccines include inactivated virus, live attenuated virus, mRNA-based, viral vectored and adjuvanted protein-subunit-based vaccines. All of them contain the spike glycoprotein as the main immunogen and result in reduced disease severity upon SARS-CoV-2 infection. While we and others have shown that mRNA-based vaccination reactivates pre-existing, cross-reactive immunity, the effect of vector vaccines in this regard is unknown. Here, we studied cellular and humoral responses in heterologous adenovirus-vector-based ChAdOx1 nCOV-19 (AZ; Vaxzeria, AstraZeneca) and mRNA-based BNT162b2 (BNT; Comirnaty, BioNTech/Pfizer) vaccination and compared it to a homologous BNT vaccination regimen. AZ primary vaccination did not lead to measurable reactivation of cross-reactive cellular and humoral immunity compared to BNT primary vaccination. Moreover, humoral immunity induced by primary vaccination with AZ displayed differences in linear spike peptide epitope coverage and a lack of anti-S2 IgG antibodies. Contrary to primary AZ vaccination, secondary vaccination with BNT reactivated pre-existing, cross-reactive immunity, comparable to homologous primary and secondary mRNA vaccination. While induced anti-S1 IgG antibody titers were higher after heterologous vaccination, induced CD4+ T cell responses were highest in homologous vaccinated. However, the overall TCR repertoire breadth was comparable between heterologous AZ-BNT-vaccinated and homologous BNT-BNT-vaccinated individuals, matching TCR repertoire breadths after SARS-CoV-2 infection, too. The reasons why AZ and BNT primary vaccination elicits different immune response patterns to essentially the same antigen, and the associated benefits and risks, need further investigation to inform vaccine and vaccination schedule development.

Host and microbiome features of secondary infections in lethal covid-19.

Secondary infections contribute significantly to covid-19 mortality but driving factors remain poorly understood. Autopsies of 20 covid-19 cases and 14 controls from the first pandemic wave complemented with microbial cultivation and RNA-seq from lung tissues enabled description of major organ pathologies and specification of secondary infections. Lethal covid-19 segregated into two main death causes with either dominant diffuse alveolar damage (DAD) or secondary pneumonias. The lung microbiome in covid-19 showed a reduced biodiversity and increased prototypical bacterial and fungal pathogens in cases of secondary pneumonias. RNA-seq distinctly mirrored death causes and stratified DAD cases into subgroups with differing cellular compositions identifying myeloid cells, macrophages and complement C1q as strong separating factors suggesting a pathophysiological link. Together with a prominent induction of inhibitory immune-checkpoints our study highlights profound alterations of the lung immunity in covid-19 wherein a reduced antimicrobial defense likely drives development of secondary infections on top of SARS-CoV-2 infection.

Hijacking of transcriptional condensates by endogenous retroviruses.

Most endogenous retroviruses (ERVs) in mammals are incapable of retrotransposition; therefore, why ERV derepression is associated with lethality during early development has been a mystery. Here, we report that rapid and selective degradation of the heterochromatin adapter protein TRIM28 triggers dissociation of transcriptional condensates from loci encoding super-enhancer (SE)-driven pluripotency genes and their association with transcribed ERV loci in murine embryonic stem cells. Knockdown of ERV RNAs or forced expression of SE-enriched transcription factors rescued condensate localization at SEs in TRIM28-degraded cells. In a biochemical reconstitution system, ERV RNA facilitated partitioning of RNA polymerase II and the Mediator coactivator into phase-separated droplets. In TRIM28 knockout mouse embryos, single-cell RNA-seq analysis revealed specific depletion of pluripotent lineages. We propose that coding and noncoding nascent RNAs, including those produced by retrotransposons, may facilitate 'hijacking' of transcriptional condensates in various developmental and disease contexts.

Generation of glycan-specific nanobodies.

The development of antibodies that target specific glycan structures on cancer cells or human pathogens poses a significant challenge due to the immense complexity of naturally occurring glycans. Automated glycan assembly enables the production of structurally homogeneous glycans in amounts that are difficult to derive from natural sources. Nanobodies (Nbs) are the smallest antigen-binding domains of heavy-chain-only antibodies (hcAbs) found in camelids. To date, the development of glycan-specific Nbs using synthetic glycans has not been reported. Here, we use defined synthetic glycans for alpaca immunization to elicit glycan-specific hcAbs, and describe the identification, isolation, and production of a Nb specific for the tumor-associated carbohydrate antigen Globo-H. The Nb binds the terminal fucose of Globo-H and recognizes synthetic Globo-H in solution and native Globo-H on breast cancer cells with high specificity. These results demonstrate the potential of our approach for generating glycan-targeting Nbs to be used in biomedical and biotechnological applications.

Sleep neuron depolarization promotes protective gene expression changes and FOXO activation.

Sleep is an essential state that allows for recuperation and survival processes. Disturbing sleep triggers stress responses that promote protective gene expression. Sleep and its deprivation grossly impact gene expression, but little is known about how normal or disturbed sleep control gene expression. Central to the induction of sleep are sleep-active neurons, which inhibit wakefulness and promote survival. Sleep and sleep-active neurons are highly conserved. In Caenorhabditis elegans, the sleep-active RIS neuron is crucial for sleep and survival. Here, we show that RIS depolarization promotes the protective gene expression response that occurs during developmental arrest. This response includes the activation of FOXO/DAF-16 and expression of DAF-16 target genes such as HSP-12.6, a small heat-shock protein that is required for starvation survival. Disturbing sleep by mechanical stimulation increases RIS depolarization. RIS activation in turn activates DAF-16 and other genes required for survival. Hence, during normal sleep, RIS depolarization promotes protective gene expression. When sleep is disturbed, protective gene expression gets further increased by raised RIS depolarization. We thus link sleep-active neuron depolarization to protective gene expression changes and suggest that the cellular stress response following sleep deprivation could be understood as a safeguarding process that is caused by the overactivation of sleep-active neurons.

Single-cell sequencing of human midbrain reveals glial activation and a Parkinson-specific neuronal state.

Idiopathic Parkinson's disease is characterized by a progressive loss of dopaminergic neurons, but the exact disease aetiology remains largely unknown. To date, Parkinson's disease research has mainly focused on nigral dopaminergic neurons, although recent studies suggest disease-related changes also in non-neuronal cells and in midbrain regions beyond the substantia nigra. While there is some evidence for glial involvement in Parkinson's disease, the molecular mechanisms remain poorly understood. The aim of this study was to characterize the contribution of all cell types of the midbrain to Parkinson's disease pathology by single-nuclei RNA sequencing and to assess the cell type-specific risk for Parkinson's disease using the latest genome-wide association study. We profiled >41 000 single-nuclei transcriptomes of post-mortem midbrain from six idiopathic Parkinson's disease patients and five age-/sex-matched controls. To validate our findings in a spatial context, we utilized immunolabelling of the same tissues. Moreover, we analysed Parkinson's disease-associated risk enrichment in genes with cell type-specific expression patterns. We discovered a neuronal cell cluster characterized by CADPS2 overexpression and low TH levels, which was exclusively present in idiopathic Parkinson's disease midbrains. Validation analyses in laser-microdissected neurons suggest that this cluster represents dysfunctional dopaminergic neurons. With regard to glial cells, we observed an increase in nigral microglia in Parkinson's disease patients. Moreover, nigral idiopathic Parkinson's disease microglia were more amoeboid, indicating an activated state. We also discovered a reduction in idiopathic Parkinson's disease oligodendrocyte numbers with the remaining cells being characterized by a stress-induced upregulation of S100B. Parkinson's disease risk variants were associated with glia- and neuron-specific gene expression patterns in idiopathic Parkinson's disease cases. Furthermore, astrocytes and microglia presented idiopathic Parkinson's disease-specific cell proliferation and dysregulation of genes related to unfolded protein response and cytokine signalling. While reactive patient astrocytes showed CD44 overexpression, idiopathic Parkinson's disease microglia revealed a pro-inflammatory trajectory characterized by elevated levels of IL1B, GPNMB and HSP90AA1. Taken together, we generated the first single-nuclei RNA sequencing dataset from the idiopathic Parkinson's disease midbrain, which highlights a disease-specific neuronal cell cluster as well as 'pan-glial' activation as a central mechanism in the pathology of the movement disorder. This finding warrants further research into inflammatory signalling and immunomodulatory treatments in Parkinson's disease.

The metabolic growth limitations of petite cells lacking the mitochondrial genome.

Eukaryotic cells can survive the loss of their mitochondrial genome, but consequently suffer from severe growth defects. 'Petite yeasts', characterized by mitochondrial genome loss, are instrumental for studying mitochondrial function and physiology. However, the molecular cause of their reduced growth rate remains an open question. Here we show that petite cells suffer from an insufficient capacity to synthesize glutamate, glutamine, leucine and arginine, negatively impacting their growth. Using a combination of molecular genetics and omics approaches, we demonstrate the evolution of fast growth overcomes these amino acid deficiencies, by alleviating a perturbation in mitochondrial iron metabolism and by restoring a defect in the mitochondrial tricarboxylic acid cycle, caused by aconitase inhibition. Our results hence explain the slow growth of mitochondrial genome-deficient cells with a partial auxotrophy in four amino acids that results from distorted iron metabolism and an inhibited tricarboxylic acid cycle.

Altered DNA Methylation Profiles in SF3B1 Mutated CLL Patients.

Mutations in splicing factor genes have a severe impact on the survival of cancer patients. Splicing factor 3b subunit 1 (SF3B1) is one of the most frequently mutated genes in chronic lymphocytic leukemia (CLL); patients carrying these mutations have a poor prognosis. Since the splicing machinery and the epigenome are closely interconnected, we investigated whether these alterations may affect the epigenomes of CLL patients. While an overall hypomethylation during CLL carcinogenesis has been observed, the interplay between the epigenetic stage of the originating B cells and SF3B1 mutations, and the subsequent effect of the mutations on methylation alterations in CLL, have not been investigated. We profiled the genome-wide DNA methylation patterns of 27 CLL patients with and without SF3B1 mutations and identified local decreases in methylation levels in SF3B1mut CLL patients at 67 genomic regions, mostly in proximity to telomeric regions. These differentially methylated regions (DMRs) were enriched in gene bodies of cancer-related signaling genes, e.g., NOTCH1, HTRA3, and BCL9L. In our study, SF3B1 mutations exclusively emerged in two out of three epigenetic stages of the originating B cells. However, not all the DMRs could be associated with the methylation programming of B cells during development, suggesting that mutations in SF3B1 cause additional epigenetic aberrations during carcinogenesis.

Dnmt1 has de novo activity targeted to transposable elements.

DNA methylation plays a critical role during development, particularly in repressing retrotransposons. The mammalian methylation landscape is dependent on the combined activities of the canonical maintenance enzyme Dnmt1 and the de novo Dnmts, 3a and 3b. Here, we demonstrate that Dnmt1 displays de novo methylation activity in vitro and in vivo with specific retrotransposon targeting. We used whole-genome bisulfite and long-read Nanopore sequencing in genetically engineered methylation-depleted mouse embryonic stem cells to provide an in-depth assessment and quantification of this activity. Utilizing additional knockout lines and molecular characterization, we show that the de novo methylation activity of Dnmt1 depends on Uhrf1, and its genomic recruitment overlaps with regions that enrich for Uhrf1, Trim28 and H3K9 trimethylation. Our data demonstrate that Dnmt1 can catalyze DNA methylation in both a de novo and maintenance context, especially at retrotransposons, where this mechanism may provide additional stability for long-term repression and epigenetic propagation throughout development.

Integrated analysis of Xist upregulation and X-chromosome inactivation with single-cell and single-allele resolution.

To ensure dosage compensation between the sexes, one randomly chosen X chromosome is silenced in each female cell in the process of X-chromosome inactivation (XCI). XCI is initiated during early development through upregulation of the long non-coding RNA Xist, which mediates chromosome-wide gene silencing. Cell differentiation, Xist upregulation and gene silencing are thought to be coupled at multiple levels to ensure inactivation of exactly one out of two X chromosomes. Here we perform an integrated analysis of all three processes through allele-specific single-cell RNA-sequencing. Specifically, we assess the onset of random XCI in differentiating mouse embryonic stem cells, and develop dedicated analysis approaches. By exploiting the inter-cellular heterogeneity of XCI onset, we identify putative Xist regulators. Moreover, we show that transient Xist upregulation from both X chromosomes results in biallelic gene silencing right before transitioning to the monoallelic state, confirming a prediction of the stochastic model of XCI. Finally, we show that genetic variation modulates the XCI process at multiple levels, providing a potential explanation for the long-known X-controlling element (Xce) effect, which leads to preferential inactivation of a specific X chromosome in inter-strain crosses. We thus draw a detailed picture of the different levels of regulation that govern the initiation of XCI. The experimental and computational strategies we have developed here will allow us to profile random XCI in more physiological contexts, including primary human cells in vivo.